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The application of linear analysis in determining a guidance func- 
tion is investigated. The differential equations of motion are linearized 
about a nominal calculus of variations solution. The result is an explicit 
expression for the cutoff radius error, Ac , and cutoff angle error, ^AB, as 
a linear operation on deviations in initial conditions and several non- 
linear functions of thrust angle deviations and thrust acceleration devia- 
tions along the trajectory. With this expression available, a suitable 
form is selected for a function to determine thrust angle, X. The coef- 
ficients of this function are mathematically determined from the explicit 
solution obtained for Ac and A0 under the constraint that these values 
be as near zero as feasible for deviations in initial conditions and 
thrust acceleration whose values are arbitrary within their expected 
range of variation. 

The results of employing this function to determine X for a number 
of examples are shown. These results emphasize the advantage of mathe- 
matically imposing the mission criteria in determination of guidance 
coefficients as well as illustrate the value of linearization techniques 
in guidance analysis. 
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TECHNICAL MEMORANDUM X-53166 


GUIDANCE APPLICATIONS OF LINEAR ANALYSIS 
SUMMARY 


An explicit solution is obtained to the linearized differential 
equations of motion. 1 To illustrate the usefulness of such a solution 
it is used to impose the mission criteria in determining coefficients for 
a guidance function. The resulting function was employed in an available 
computer program designed to determine guidance function performance. 

The results were exceptionally good and serve to emphasize the advantages 
of this type of analysis. 


INTRODUCTION 


To determine a forcing function in such a way that it accomplishes a 
given result, it is useful to first determine the effect that this forcing 
function and other parameters and forcing functions have on this result. 
With this knowledge available, the task is considerably simplified. This 
report applies this principle to the problem of guidance. The equations 
of motion are linearized about a nominal trajectory. These equations are 
solved so that the important cutoff deviations, Ac and A0, are obtained 
explicitly as a function of initial conditions, thrust acceleration and 
thrust angle, A form is selected for the function which determines the 
thrust angle, X. The explicit solution for At and A0 is employed to 
impose the mission criteria in determining the coefficients of this 
function which, as expected, accomplishes well the result for which it 
was designed. The method employed to obtain the explicit solution, the 
manner in which it was employed to determine the guidance function and the 
results obtained from using this guidance function in a number of examples 
are presented and discussed in detail. 

The explicit solution to the linearized differential equations of 
motion was programmed on the IBM 1620 computer by Mr. Quintin Peasley, 
Technical and Scientific Staff, Aero-As trodynamics Laboratory, George C. 
Marshall Space Flight Center. 


1 A similar analysis was effectively employed in Reference 3. The 
accuracy of the solution and the results it produced led directly to 
conclusions that would not likely have otherwise been suspected. 



CHAPTER I. EXPLICIT SOLUTION 
Section 1 # Linearized Equations of Motion 

The explicit solution desired requires the solution to the linearized 
equations which describe the motion of the vehicle. For this analysis, 
the motion was assumed to be in a plane and described by the following 
differential equations : 


x = - s in X + x . 
m g 


y = — cos X + y 
y m 2 


( 1 . 1 . 1 ) 

( 1 . 1 . 2 ) 


The solution to the above set of differential equations is assumed 
to be known numerically for a standard set of conditions. The question 
concerning the solution under nonstandard conditions is now considered. 
Employing no approximations yet, the equations under nonstandard condi- 
tions can be written as a function of the standard solution and deviations 
from the standard as follows: 


x + Ax = (- + — Vs in X + A Sin x) + x + Ax . 

Vm m r g g 


y + Ay = + — j (cos X + Acos X) + y g + Ay g . 


These can be rearranged to form 


x + Ax = | Sin X + X g + ^1 + 77^) | 


A sin X + sin X — + ^X 
m e 


y + Ay = ~ cos X + y + 
J J m g 


i + 


Af/m 

f/m 


— A cos X + cos X — + Ay 

m mg 


3 



Substitution of equations (1.1.1) and (1.1.2) into the above expres- 
sions yields the following exact expressions: 


Ax = + ( 1 + ^7“^) A sin X + sin X — 


f /m 


/ 


m 


Ay = Zy„ + ( 1 + A cos X + cos X 


At this point, the first approximation is employed to determine 
£x'g and Xyg as functions of position deviations Zsc and Zy. 


Zx = h, ZiJC + h 2 Zsy. 
g 1 


(1.1.3) 



= k x ZJi + k 2 zy, 


where 


and 



(1.1.4) 


The next approximations can be extended to include as many terms 
as necessary. The three terms included below were found sufficient for 
this report. 


A s in X = 


TC 

180 


cos X AX 





4 





With the following definitions, the system can be expressed in 
more convenient matrix notation. 



5 







In this notation, the differential equations can be written 


AX = AAX + HAF, 


(1.1.5) 


where 


AX = 


Ax 

Ay 

Ax 

A7 


The terms Ac and Zy are expressed in km, and /^y in m/sec, AX in 
degrees, and Af/m in m/sec 2 . 

The time variable deviations required in equation (1.1.5) are 
defined as 


AX = X(t + At 0 ) - X s (t) 
AX = X(t + At 0 ) - Xg(t) 
Af/m = | (t + At 0 ) - |(t) 
At 0 = tf - to. 


where tj. is second stage ignition time on any trajectory, t Q is second 
stage ignition time on the standard trajectory, and the subscript s 
refers to values obtained from the standard trajectory. 
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Section 2. Explicit Solution to the Differential Equations 


The solution to equation (1.1.5) is shown in Reference 1 to be of 
the following form: 


AX n = U(t n , 


to) 


U(t n , t) H(t) AF (t) dt, 


( 1 . 2 . 1 ) 


where U(t n , t) is the solution to the differential equation below 
evaluated at t x = t n from initial conditions at t. 


U(t lf t) = A(t x ) U(t l3 t), U (t , t) = I. (1.2.2) 


The solution to equation (1.2.2) is obtained by assuming that the 
elements of the A matrix are constant over each of a number of small 
intervals. This assumption yields the solution at the end of this small 
interval which is then used to provide initial conditions from which the 
solution at the end of the next small interval is obtained. Continuing 
in this manner, the solution U(t n , t) can be obtained for t = t Q and any 
of a number of values of t, t Q ^ t g t n , where t n is cutoff time on the 
s tandard trajectory. 

Under the assumption just made, the differential equations which 
will be solved are 


AX(t) = A(£ k ) ZSC(t) + H(t) AF( t), t k _ x ^ t g t k 


where 




\ + Vi 


k = 1, 2, 


n. 
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For this system, the solution to equation (1.2.2) can be written as 


U(t, 


fc k- ^ 


= A(| k )At k 


A (tJ At i 




i=0 


(1.2.3) 


This solution appears in Reference 1 and is adapted to general 
application to large systems in Reference 2. Truncation error associated 
with the series in equation (1,2.3) is sometimes a problem. However, 
the A matrix considered in this application has such small elements and 
is changing so slowly with time that At^ ^ 5 sec was found adequate and 
only the following terms were included in the series approximation: 


e A(! k )At k = j + A( | k ) A t k . 


Then 


U(t 


n’ 


t 

n-i 


) = 


e A (ln) At n 


and 


U(t n> ’Hc-i 5 = U(t n’ H) e A(lk)Atk > k = n - 1, ..., 1. 


In this manner, U(t n , t 0 ) can be evaluated numerically and U(t n , t) 
for any value of t desired. Although the integration indicated in equa- 
tion (1,2.1) might require more sophisticated techniques, the functions 
encountered in this application are sufficiently smooth that this inte- 
gral can be evaluated quite well by a summation, where each element of 
the sum consists of the integrand evaluated at the midpoint of an 
interval t k _ x g t g t k multiplied by the length of the interval At k for 
values of At k as large as forty seconds. For extreme accuracy, however, 
At k was chosen at values on the order of 5 seconds. 
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Section 3, Extrapolation to Cutoff Time 


The procedure just described provides the solution at standard 
cutoff time t n . The following equation will be used to determine AX C , 
the deviations in state variables obtained by subtracting their standard 
value at t n from their actual value at a different cutoff time, t c : 


Z^X = AX + 
c n 



x(t + At 0 ) dt, 


( 1 . 3 . 1 ) 


where X(t 4- At Q ) is evaluated on the nonstandard trajectory. This can 
be written 

X(t + At 0 ) = X s (t) + X(t + At 0 ) - X s (t) 

X(t + At ) = X (t) + AX(t). 

O o 

These terms can be further decomposed: 


t 

X (t) = X + / X (t) dt. 

s n J s 

t 

n 


£X(t) - AZX(t) + H(t) AF(t). 


Subscript s refers to values obtained from the standard trajectory. 
The subscript n refers to values at cutoff time, t n , on the standard 
trajectory. 
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can be written as 


The term AX C 


t -At 
c o 


AX = AX + X At + 
c n n 


I dt * 


n 



X (t) dt + 
s 



[AAX 4- HAF] dt, 


where 


At = t 

e 


t 

n 


At . 
o 


This can be further simplified, by notation, to the following expression: 


At = X At + E, 
c n 


(1.3.2) 


where 


t -At 
c o 


t -At 
c o 


E = At 


n + J dt x J x g (t) dt + J [AAX + HAF] dt. (1.3.3) 


n 


The value of At in equation (1.3.2) depends on the cutoff criterion. 
Under the assumption of cutoff at a constant velocity, the following 
will be used to determine At. 


AV 

c 




= 0 , 


or 


AV = T, AX„ = 0, 
c c 


(1.3.4) 


where 



n 


0 x 

n 



(1.3.5) 
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Equations (1.3.2) and (1.3.4) give the following relationship: 


Tn AX = T, X At + T-iE = 0, 
x c 1 n x 


and 


At 


T X E 


T, X 


(1.3.6) 


Substitution of this expression for At into equation (1.3.2) yields 


AX 

c 


X T n E 


T x X 


n 


+ E, 


or 


AX = T E, (1.3.7) 

c o 

where 

X T x 

T = I - -2-_ , (1.3.8) 

O n, ,, 


and E is defined in equation (1.3,3). This expression for E requires 
some assumptions or approximations in order to actually be evaluated. 
The following approximations are used: 


• . • • 

X (t) = X 
s n 

and 


AAX(t) + H(t) = A Z& + H £F 
nn nn 
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for 



3 

IIA 

t s t n 

+ 

At, 

At 

> 

+ 

n 

fc 

flA 

rt 

< 

t n’ 

At 

< 


o- 


0 - 


These approximations essentially keep second order terms in the 
expression for E, discarding only terms of third order or higher. The 
resulting expression for E with these approximations is 


E = At + X 
n n 


At* 


+ A AX At + H At AF . 
n n n n 


(1.3.9) 


Equations (1.3.7), (1.3,8), and (1,3.9) give At and AX C as functions 
of AX n . Equation (1,2,1) yields the solution of AX n as a function of 
initial conditions, At 0 , deviations in thrust acceleration, Af/m, and 
deviations in thrust angle, AX, for toStS t n . The other functions 
appearing in these equations are evaluated from data on the standard 
trajectory. Thus, with a given standard trajectory, an explicit expres- 
sion for cutoff deviations for nonstandard trajectories has been obtained. 


Section 4, End Conditions 


To meet the mission, certain end conditions are required. The 
particular variables for which an explicit solution is to be obtained 
depend on the mission. Under the assumption of an orbital mission, Ar 
and A9 are the variables of concern. The following approximations can 
be used to determine these variables as a function of the vector At c : 

x y 

Ar = — As + — A7 (1.4.1) 

r . c r c 
n n 


A < r c v c cos e c ) = A(x c x c + y c y c ). 


A(r v ) cos e + r v A cos 9=Axx+xAx+^yy+yZ^y. 
cc n nn c cn n c c n n c 
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Although the following assumption is not necessary, it is conven- 
ient at this point to take advantage of the fact that the mission under 
consideration is a circular orbit at a fixed radius. Then r c is constant 
and A(r c V c ) = 0. Then the following approximation is used: 


A cos 8 = 

c 


180 si " e n Ae c' 


Since sin 0 n = 1, we have 


A cos 0 = 

c 


180 A9 c' 


The following expression is then obtained for A© c * 


180 

m v” < A „ + + x n + y n ^ c> ‘ 

n n 


(1.4.2) 


Equa tions 
equation. 


(1.4,1) and (1.4.2) can then be combined into one matrix 


AR ■" Tp , 

^ c 


(1.4.3) 


where 




Ar 

Ae c 


AG 

-J 


— 


and 



-180* 

n 



v 

n 



r 

n 


-180y n 


VX 


n 


v 

n 


-180x 


n 


vX v 
n n 


- I80y 


jtr v 
n n 


(1.4.4) 
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Equations (1.2.1), (1.3.7), (1.3.8), (1.3.9), (1.4.3) and (1.4.4) 
give us the means of determining the variables At and A0 as a function 
of initial condition deviations, Z^q, thrust acceleration deviations, 

“(t), and thrust angle deviations, £X(t), (t Q g t g t n ) . This expres- 
sion can be used to actually evaluate errors of individual trajectories 
in the neighborhood of the standard if all the deviations mentioned are 
known. It also provides considerable insight into the mechanics of 
solving differential equations of motion by showing, term by term, the 
effect on mission error of 

AK q , •~( t )> and AX.(t). 


In addition, for this particular problem, the explicit representation 

provides a means of determining X as a function of AX 0 and “(t) so that 
Ac and A9 are as near zero as this analysis and the assumptions concern- 
ing the form of X will allow. Although AX OJ “—(t), anc * AX(t) were the 

only parameters considered in this analysis, with little additional effort 
other forcing functions or parameters could have been included. For the 
present, the concern will remain with the forcing functions and parameters 
already considered and an actual application will be demonstrated. 


Section 5. Numerical Example 

Several calculus of variations solutions were available for an 
early SA-6 second stage vehicle. The standard trajectory had the follow- 
ing initial and end point conditions: 


Initial Conditions : (1.5.1) 


X = 153.98343 

km 

t = 146,815 sec 

o 


o 

y D = 6435.8783 

km 

f _ 8.78065 

m U; 1.3751 - . 20888x 5 

x Q = 2818.3294 

m/sec 

t - t 



_ o 

y Q = 988.35767 

m/sec 

T 100 • 
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Final Conditions: 


(1.5.2) 


x = 2326.37 km 
n 


v = 7792 m/sec 
n 


y = 6128.51 km 
n 


r = 6555.200 km 
n 


x = 7285.34 m/sec 
n 


t = 620.679 sec 
n 


y = - 2765.50 m/sec 
n 


x = 19.047 m/sec 2 
n 


y = - 13.596 m/sec 2 
n 


*x - .1001 m/sec 3 
n 


y* = - .0655 m/sec 3 


-(t n ) = 22.793 m/sec‘ 
m 11 


X = 102.506 c 
n 


e n = 90°. 


The gravity components were defined by 


•• x ~ 

x = — g 
g r 


= Z 


where 


g = - 


2 

g o r 0 


g n - 9.81 m/sec 2 , 


r = 6370 km. 
o 


From this, the following elements of the A matrix are determined. 
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[1 - 3(x/r) 2 ] 


hi 



5 r 2 
o o 


1 


h 2 




ki 



g r 2 
o o 



k 2 



g r 2 
o o 


[1 


3(y/r) 2 ] 

/ 


(1.5.3) 


The values of x, y, and t obtained from the standard trajectory are 
listed in Table 1.1. They are listed as a function of t where t desig- 
nates time on the standard trajectory. In addition, the quantity At^ is 
listed which describes the length of the interval over which the dif- 
ferential equations were assumed to be a constant coefficient system. 


TABLE 1.1 

STANDARD TRAJECTORY DATA 


t(sec) 

T(10 2 sec) 

Atk(sec) 

x(km) 

y (km) 

150 

.0318 

13.18 

163.0 

6435.6 

180 

.3318 

40 

250.4 

6461.9 

220 

.7318 

40 

374.8 

6488.2 

260 

1.1318 

40 

508.6 

6504.4 

300 

1.5318 

40 

652.7 

6510.5 

340 

1.9318 

40 

807.8 

6506.2 

380 

2.3318 

40 

975.1 

6491.0 

420 

2.7318 

40 

1155.7 

6464.5 

460 

3.1318 

40 

1351.0 

6426.0 

500 

3.5318 

40 

1562.7 

6374.5 

540 

3.9318 

40 

1793.0 

6308.8 

580 

4.3318 

40 

2044.7 

6227.1 

610 

4.6318 

20.68 

2249.9 

6153,9 
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The data in Table 1.1 can be used to evaluate the elements of the 
A matrix defined in equations (1,5.3). The results are shown in 
Table (1.2) which for convenience have been multiplied by 1000, 

TABLE 1.2 


ELEMENTS OF A MATRIX 


t 

h i 

h P 

k n 

k P 

150 

-1.489 

.113 

.113 

2.327 

180 

-1.465 

.171 

.171 

2.937 

220 

-1.436 

.251 

.251 

2.886 

260 

-1.407 

.334 

.334 

2.840 

300 

-1.378 

.423 

.423 

2.799 

340 

-1,348 

.518 

.518 

2.761 

380 

-1.315 

.621 

.621 

2.722 

420 

-1.275 

.731 

.731 

2.682 

460 

-1.228 

.849 

.849 

2.633 

500 

-1.168 

.977 

.977 

2.576 

540 

-1.113 

1.079 

1.079 

2.505 

580 

-1.001 

1.257 

1.257 

2.415 

610 

- .912 

1.370 

1.370 

2.327 


NOTE: All of the above elements have been multiplied by 1000. To use 

in the construction of the A matrix, they must first be multiplied 
by 10“ 3 . 


The matrix 
Tables 1.1 and 


U(t n , t) can be determined from the information in 
1.2 for the values of t listed as follows. 


u(t r 


t ) = I + A(! ) At 
n~i' v n' n 
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where 


0 


0 10" 3 0 


A U n > - 


Ml n > 

k i^n> 


h ^n> 

M£ n ) 


o 10" 3 

o o 

o o 


Choosing | n = 610 and t n = 600, t n = 620,68, then the following matrix 
is determined. 


U(t n , 600) = 


1 

0 

,01886 

,02833 


0 

1 

.02833 

.04812 


.02068 

0 

1 

0 


0 

.020681 
0 
1 


Then U(t n , 560) = U(t n , 600) U(600, 560) where U(600, 560) = I + 40A(580). 

Continuing in this manner, U(t n , t) can be determined for a number 
of values of t back to and including t = t Q . The following matrix 
U(t n , t Q ) was obtained from the data in Tables 1.1 and 1.2. 


.86115 

.04792 

.45637 

.00905 

.04503 

1.29901 

.00890 

.51145 

.57740 

.35385 

.88097 

.09625 

.31006 

1.38267 

.09284 

1.27619 


(1.5.4) 


The intermediate matrices U(t n , t) are tabulated in Appendix I. 
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The expression for U(t^, tk- 1 ) is the truncation of a series expan- 
sion. If the elements in the A matrix are large, more terms are required. 
If they are extremely large, U(tk, t^-.. j) must be evaluated by special 
methods described in Reference 2. In addition, the interval over which 
the system is assumed to have a constant A matrix depends on the particular 
problem at hand. The numerical example being used as an illustration 
assumes a constant coefficient system for intervals of forty seconds. 

A more accurate solution was employed to determine the guidance function 
in Chapter II. The system was assumed constant for intervals of 5 seconds 
or less, and the integration was carried out by summing over five second 
intervals. The solution differed generally in the third significant 
figure from that obtained from the present example where 40 seconds was 
used . 


In addition to the matrix U(t n , t Q ), the solution given by equation 
(1.2.1) requires that an integration be performed. The functions in this 
example are sufficiently smooth that a very good solution can be obtained 
by summing U(t n , tk) H(t^.) AFCt^) At^.. If the functions in the integrand 
are extremely variable, this technique may not give accurate values for 
the integral. Although techniques outlined in Reference 2 overcome this 
problem, they are not necessary for this example. With the following 
definition, AX n can be written as a linear sum of deviations in initial 
conditions AX 0 and the vector AF(t) evaluated for several values of t. 


U (t . ) = U(t , 
J n’ 


t.) H ( t ; ) At., 
J J J 


(1.5.5) 


Then, 


AX = U(t , 
n n 


t ) AX 
o o 


U(t ) AF(tj) . 


(1.5.6) 


j = l 


Table 1.3 lists the data taken from the standard trajectory 
necessary to evaluate the matrix H(t*) for the values of tj indicated. 

In addition, the value of Atj is listed which was used to determine 
U(tj) as defined in equa tion J (1 . 5 . 5 ) . Appendix I lists the corresponding 
matrices U(t n , tj) which, with the data in Table 1.3, is sufficient to 
evaluate the elements of U(tj) necessary to perform the summation indi- 
cated in equation (1.5.6). 
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TABLE 1.3 


DATA REQUIRED FOR H(tj) Atj 


t . (sec) 
J 

At . (sec) 
1 

X (°) 

— (m/sec 2 ) 
m 

160 

33.18 

59.255 

6.516 

200 

40 

62.460 

6.946 

240 

40 

65.745 

7.438 

280 

40 

69.119 

8.005 

320 

40 

72.590 

8.665 

360 

40 

76.168 

9.444 

400 

40 

79.861 

10.376 

440 

40 

83.676 

11.513 

480 

40 

87.620 

12.930 

520 

40 

91.693 

14.743 

560 

40 

95.896 

17.149 

600 

40.68 

100.223 

20.493 


The above values, together with U(t n , t i ) listed in Appendix I, 
were used to determine U(tj) from equation (1,5 # 5). The elements of 
U ( 1 4 ) are listed in Appendix II. All of the information necessary to 
evaluate ZiXn available and the necessary coefficients evaluated; 

AX 0 , AX(tj) and —(tj) remain explicit. 

The last step necessary to obtain an explicit expression for Ar and 
A0 is to evaluate the matrices defined in equations (1.3.7), (1.3.9), and 
(1.4.3). From the end' conditions given in (1.5.2), T x , defined in equa- 
tion (1.3.5), is evaluated. 

T x = (0 0 .9350 -.3549). (1.5.7) 
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Evaluation of the elements in T 2 , defined by equation (1.4.4), 
gives the following result. 


T 2 = 


.3549 .9349 0 

.008172 .003102 -.002609 


-.006875 



7.2853 


.01905 


-2.7655 

and X = 

-.01360 


19.047 

n 

.1001 


-13.596 


-.0655 


Equation (1.3,7) provides the relationship, 


AX = T E, 
c o 


where 



Equation (1.4.3) gives the relationship 

AR = Tp AX • 

^ c 

Combining these expressions yields 

/ * T i\ 

AR = T 2 ( I - — ) E = TE, 

’ v T 1 V 


(1.5.8) 
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where 


T = 



X T- 


T-, X 


(1.5.9) 


Numerical evaluation of the elements gives the following matrix T. 


3549 

.9349 

0 

0 






(1.5.10) 

008172 

.003102 

-.001603 

-.007262 



The final expression for E is given by equation (1.3.9). 


E = AX + X 


At* 


n 


+ A ZNX At 
n n 


H At AF . 
n n 


In addition to the quantities already evaluated, A n and H n must be 
evaluated in order to investigate the second order terms. From cutoff 
data, given by (1.5.2), the following can be determined. 


A 


n 


0 

0 

-.882 

1.407 


0 

0 

1.407 

2.297 


1 

0 

0 

0 


0 

1 

0 

0 


x lO' 3 


H = 
n 


0 

0 

-.0853 

-.3870 


0 

0 

0 

0 

0 

0 

0068 

.0000 

-.2152 

0015 

.0001 

.9766 


The effect of these matrices on AR is determined by equation (1.5.8) and 
(1.3.9). 
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m = TE = TAX n + 


TX —■ + TA Z& + TH £F 
n Z n n n n 


At. 


(1.5.11) 


Equation (1,3.6) gives the following expression for At, 


TtE 

At - - “ 1 — 


T X 


Using the first order approximation for E, this expression becomes 


At 


T 1 

^ n 


Ti 


X 


n 


This can be substituted for At inside the brackets in equation (1.5.11), 
and the following expression results: 


AR * T£X + 
n 


T A - 
n 


X n 

2T l X n 


AX + TH AE 
n n n 


At. 


(1.5.12) 


A brief look at the values of the elements of the matrices inside the 
brackets will give an estimate of the contribution of these second order 
terms . 

Cutoff data from the standard trajectory, equations (1.5.2), pro- 
vide the information necessary to evaluate these matrices. 



X T, 


2Tj. X r 


0 


-.0088 


0 


-.0189 


.4779 


-.0106 


.8882 

.0040 


x 10" 3 . 
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To get an idea of the effect this term might have some extreme values 
can be assumed for AX^ First, it will be assumed that Ax n - Ay n = 100 km 
and = - A£n = 100 m/sec. Then, 



X T n 
n i 


2Ti X n ' 




/-4l , 02\ 

( ) x 1CT 3 . 

\ -4.23/ 


This means that the contribution to Ar is -.041 km or -41 m for each 
second of additional second stage burning time At. However, for A9 the 
contribution is only -.004° for each second At. The signs associated 
with the deviations assumed for AX n were chosen to have the greatest 
effect on A0. With this in mind, it can be concluded that these second 
order terms are negligible with respect to AG. An extreme example for 
Ar might be AK n = Ay n = 100 km and A£ n = Ay n = 100 m/sec. This gives 



X T, 
n 1 

2 Tl X 


AX„ 


A36.6K 

( ) x 10" 3 . 

X -3.43/ 


This gives a smaller error in A0, but in Ac it would produce about 137 m 
for each second of additional second stage burning time At. Whether this 
term is negligible or not depends on a more critical investigation of 
Ac n and A/n which are not independent of each other. In most of the 
cases considered later, Z& n and Z$ n were of opposite sign so that their 
effect on Ar tended to cancel rather than add. A more detailed examina- 
tion of the relationship of Ax n and Ay n and their effect on the error in 
Ar will not be pursued in this report. However, in the event that any 
decisions depended on the assumption that this contribution be negligible, 
further investigations could be completed at that time. 

The other term in the brackets in equation (1.5.12) is 


THn 


0 

.00295 


0 0 0 

m 

.0000 .0000 -.0062 
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Choosing AX n =3° and 


gives 


Then 


“(t n > = 
m 11 


,01 -(t ) = - . 22 m/sec 2 
m n 


AF = 
^ n 


3 

9 

27 

-•22 


TH AF = 
n n 


.00885 + .0014 


The contribution from AXn is .00885° per second and from Af/m, only 
.0014° per second. The contribution of Af/m will be neglected, but the 
contribution of will be kept since it may be several hundredths of 
a degree. 

Having investigated the possible contribution of various second 
order terms and finding several of them to be insignificant, we will 
discard these and rewrite the expression for AR. 

Defining two row vectors V x and V 2 by the following equation, 


V = 




X T r 
n 1 

T i X 


then. 


AR = TAX 4- VAX At + TH AF . 

n n n n 


(1.5.13) 
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Equation (1.2,1) gave the following expression for 1 £X n , 


n 


AX n = U(t n , t Q ) + / U(t n , t) H(t) AF(t) dt. 


Substitution of this expression for £X n into the linear term T/S n appear- 
ing in (1,5.13) gives 


= TU(t , 


t 

n 

V + / TO( V 

t 

o 


t) H(t) AF(t) dt + Vz^X n At + TH^ AF^. 


For convenience, these definitions will be made: 


TU (t , 
n 



TU(t , t) H(t) 



f 2 f 3 
S 2 §3 




Since all but one of the elements in the matrix TH n were zero or 
considered negligible, the following simplification is used: 


TH AF 
n n 



AX . 

n 
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With these definitions and simplifications, the following expres- 
sions for Ar and A0 are obtained: 


= Un AX + F AF dt + V, AX At 

- 1 o / x n 


A0 = U 2 AX + / G AF dt + g 5 AX At. 

o I 3 n 


The numerical values obtained for these expressions are listed below. 


U-l = .34772 1.23145 .17029 .48137 


Uo = -.008224 -.006970 -.005788 -.007909 


V x = o 0 .4779 .8882 x 10 


g 5 = .00295. 

The term AX n is defined by equation (1.5.6). U(t n , t Q ) is found 
in Appendix I and the elements of U(tj) are listed in Appendix II. 

The first order approximation for At is 


Ti AX 


At = - 


= (0 0 -.04131 .01568) AX n , 


where AX n can be evaluated from equation (1.5.6) as already mentioned. 
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The integrals appearing in equations (1.5.14) and (1,5.15) can be 
evaluated by summation. 


n 


G AF dt = ) G (t . ) AF(t.) At. 

i_ i <3 J J 

j = l 


and 


n 


F AF dt = ) F(t.) AF(t.) At., 

(_ i J J J 

j=l 


The multiplying factors for AF(tj) will be defined as follows: 


G(t^) = G(0 Atj = (g i( fc j) g3(t j ) §4(t j ) 


F(tj) 


F(t.) At. = 
J J 








These quantities are listed in Tables 1,4 and 1,5. The Atj employed 
to determine these values is listed in Table 1.3. 
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TABLE 1.4 


F<V 


tj 

fi 

f 2 (10“ 2 ) 

f 3 (10" 4 ) 

h 

160 

-1.1877 

-1.284 

.6030 

12.930 

200 

-1.4612 

-1.388 

.7419 

13.118 

240 

-1.4790 

-1.262 

.7508 

11.144 

280 

-1.4843 

-1.133 

.7535 

9.297 

320 

-1.4741 

-1.000 

.7484 

7.581 

360 

-1.4444 

- .863 

.7333 

6.004 

400 

-1.3897 

- .723 

.7055 

4.576 

440 

-1.3023 

- .581 

.6612 

3.313 

480 

-1.1715 

- .431 

.5948 

2.219 

520 

- .9802 

- .297 

.4977 

1.319 

560 

- .7023 

- .163 

. 3566 

.624 

600 

- .2958 

- .047 

.1502 

.151 

600* 

- .2363 

- .037 

.1199 

.151 


"The values in this row were obtained by using Atj = 32.5 sec instead 
of 40.68 sec. If AX = 0 for t > 612.5, these values should be used 
instead of those listed in the regular table for tj = 600 seconds. 


30 



TABLE 1.5 


6(t ) 


t 

1 

ii ( io - 2 ) 

g 2 (10“ 4 ) 

g 3 ( io - 6 ) 

OQ 1 

160 

1.438 

2.972 

-.730 

-.2995 

200 

2 . 068 . 

3.487 

- 1.049 

-.3297 

240 

2.440 

3.436 

- 1.239 

-.3033 

280 

2.860 

3.366 

- 1.452 

-.2761 

320 

3.338 

3.271 

- 1.695 

-.2479 

360 

3.887 

3.144 

- 1.973 

-.2186 

400 

4.525 

2.973 

-2.297 

-.1881 

440 

5.276 

2.747 

- 2.679 

-.1567 

480 

6.179 

2.436 

- 3.137 

-.1237 

520 

7.289 

2.015 

- 3.701 

-.0897 

560 

8.702 

1.427 

- 4.418 

-.0547 

600 

10.764 

.582 

- 5.465 

-.0186 

600 * 

8.599 

.465 

- 4.366 

-.0186 


*The values in this row were obtained by using At. = 32.5 sec instead 
of 40.68 seconds. If AX = 0 for t > 612.5, these J values should be used 
instead of those listed in the regular table for tj = 600 seconds. 
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The results just described give the influence coefficients for AX, 
Af/m and AX 0 a s they affect AT and A0. In that respect they are of 
interest in themselves for the insight they provide. For computational 
convenience, large integration steps were taken. This was done on a 
desk calculator, and the procedure outlined can be checked against these 
results by a desk calculator to verify the steps involved. Although 
these results are accurate through two significant figures, this is not 
sufficient for later work where the assumption of a constant coefficient 
system and the integration steps were necessarily reduced from 40 second 
intervals to 5 second intervals. These more accurate results were 
obtained from an IBM 1620 digital computer program and used to derive 
the guidance function discussed in the following chapter. Nevertheless, 
the numerical example with stepwise evaluation included in this chapter 
should serve to illustrate the procedure described. 
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CHAPTER II 


DETERMINATION OF A GUIDANCE FUNCTION 


Section 1, Fitting the Nominal Trajectory 

Equations (1.5.14) and (1.5.15) of the preceding chapter provide an 

explicit solution for AT and A0 as a function of Z$: 0 , ““(t) ancl AX(t). 

These equations can be employed to ensure that X, a function derived to 
approximate the function X on the standard trajectory, will meet the 
required end conditions, Ar = A9 = 0. Under standard conditions AX Q = 0 

and ~*(t) = 0. In addition to this, if X sufficiently well approximates 

the nominal value of X, the contribution of higher order terms becomes 
negligible. It will be assumed that this can be done sufficiently well 
by describing X as a quadratic in t. For convenience, the following 
transformation will be employed. 



where t* denotes time on any trajectory and t-^ denotes ignition time on 
that same trajectory. The term t 1 on an arbitrary trajectory will be 
related to t on the standard trajectory by the following equation, 


t' = t + At os 


( 2 . 1 . 2 ) 


where 


At o = 


t. - t . 
i o 


(2.1.3) 


Substituting these expressions into equation (2.1.1) gives 


t - t 

o 


100 
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This definition was already used in equation (1.5.1). Thus, equation 
(2.1.1) offers no contradiction. 

The term X will now be defined as 


X = C o + CjT + C 2 t 2 . (2.1.4) 

With this definition and otherwise standard conditions, equations (1.5.14) 
and (1.5.15) become, respectively, 


xr 



fi 


C o + CjT + C 2 t 2 


t 

O 


x s (t) 


dt 


(2.1.5) 


A6 



C Q + C X T 4- c 2 x 2 


x s (t) 


dt. 


( 2 . 1 . 6 ) 


The higher order terms in these expressions have been considered 
negligible. The term Xg (t) is the nominal calculus of variations solu- 
tion for X. To ensure that the end conditions are met, the following 
constraints are imposed on the coefficients C 0 , and C 2 . 


C 

o 


fl dt 


+ Cq 


;n 


fj/r dt = 



f x T 2 dt 4- 



f 1 X s (t) dt. 


(2.1.7) 



r n 

r n 

r n 

r n 

~ C o 

J g i d t 4- C X 

J g x T dt = ■ C 2 

J giT 2 dt + 

J §i^ s ( t ) dt. 


( 2 . 1 . 8 ) 
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The integrals in equations (2.1.7) and (2.1.8) can be evaluated 
quite well from the coefficients in Tables 1.4 and 1.5 together with 
the values of t and Xg tabulated in Table 2.1. Although these r.esults 
may differ in the third or fourth significant figure from the values 
shown below, the values presented were those actually used. Since the 
resulting function X fits sufficiently well, the results of the following 
constraints were not corrected. 


14.284 C Q + 28.8545 C x = - 79.49615 C 2 + 1080.9536 

.5872 C + 1.80745 C, = - 6.491105 C 2 + 50.32710. 
o ± 

Solving for C D and C x yields 

C = 4.914469 C 2 + 56.523779 (2.1.9) 

o 

C ± = - 5.187907 C 2 + 9.480948. (2.1.10) 

Substituting these values into the expression for X shown in equa- 
tion (2.1.4) yields 

X = a(x)C 2 + b(x), (2.1.11) 


where 


a ( t) = 4.914469 - 5.187907x + t 2 


and 


b ( t) = 56.523779 + 9.480948x. 


Any choice of C 2 in equation (2.1.11) will produce a Ac and AG of zero 
providing the resultant value of AX. is sufficiently small that the higher 
order terms are negligible. A value of C 2 must be chosen with this in 
mind. The sum of the squared residuals are given in the following 
expression: 
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X(t.) - X(t ) . 

j=l- J 

Minimizing the value with respect to the choice of C 2 gives the following 
least squares constraint. 



n 


= 2 
Sc 2 fa 


x(t.) - x s (t.) 


1 SX(t.) 


JL-- 0. 




Since 


dX(r,) 

1_ _ a ( t . ) , 

SCo J 


this becomes 


n 


2 n 

rJT — 

Cs z 

a(T j > . 

- 1 a< V 
1 

_b( T .) - x s ( t .)_ 


j=i 


j=i 


(2.1.12) 


Using the expressions for a(-r) and b(x) indicated in equation (2.1.11) 
and the values of t • and Xg(t.) indicated in Table 2.1 gives the follow- 
ing numerical expression for equation (2.1.12): 


40.78235 C 2 = 14.57552. 


This result, together with equations (2.1.9) and (2.1.10), yields the 
following coefficients for X. 


C 0 - 58.2802 N 
Ci = 7.6268 V , 

C 2 = .3574 

I 


(2.1.13) 
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and 


X = 58.2802 + 7.6268 t + ,3574t 2 . 


(2.1.14) 


X in the above expression is given in degrees and t is defined by 
equation (2.1.1). 

Table 2.1 shows a comparison of X with Xg, the function to which 
it was fitted. The resulting residuals are sufficiently small that their 
effect on second order terms is negligible. 


TABLE '2.1 

COMPARISON OF X AND X s 


: . (sec) 
J 

(10 2 sec) 

X(t.) (°) 

vv <°> 

(X - X ) (°) 

b 

160 

.1319 

59.292 

59.225* 

.067 

200 

.5319 

62.438 

62.460 

-.022 

240 

.9319 

65.698 

65.745 

-.047 

280 

1.3319 

69.072 

69.119 

r^« 

o 

• 

1 

320 

1.7319 

72.561 

72.590 

-.029 

360 

2.1319 

76.164 

76.168 

-.004 

400 

2.5319 

79.882 

79.861 

.021 

440 

2.9319 

83.713 

83.676 

.037 

480 

3.3319 

87.660 

87.620 

.040 

520 

3.7319 

91.720 

91.693 

.027 

560 

4.1319 

95.895 

95.896 

-.001 

600 

4.5319 

100.184 

100.223 

-.039 


* The value 59.225 was erroneously used in the fit shown here. The 
actual value, as shown in Table 1.3, was 59.255. This would make the 
actual residual at this point .040 instead of .067. Since the coef- 
ficients in equation (2.1.13) were used in later work and since the 
error was of minor significance, for consistency and economy of effort, 
this correction was not made and the coefficients in X have been employed . 
as they appear in equation (2.1.14). 
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Section 2. Initial Conditions 

If we take another look at equations (1.5.14) and (1.5.15), they 
appear as follows: 




To keep the guidance function as simple as possible, the second 
order terms must be negligible or made so. The term V x AX n At is 
extremely difficult to evaluate at second stage ignition. Both At and 
AX n are functions of Af/m and AX along the entire trajectory. Because 
of the difficulties involved, this term will first be ignored, and if 
after making this assumption, the resulting function is not sufficiently 
accurate, the problem can again be taken up at that time. 
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The second order term g 5 £X. n At, however, poses a different problem. 
The forcing function which we wish to determine, AX, appears as a multi- 
plier. This term was the only one considered significant in the follow- 
ing term from the expression for E, equation (1,3.3), after it was 
substituted into equation (1.5.8). 


-At, 


-At r 


TH AF dt 


TH-j 1 + 


Af/m 

f /m 


AX dt. 


n 


n 


( 2 . 2 . 1 ) 


If AX is chosen to be zero throughout the interval defined by the limits 
of integration, the entire integral would be zero. This can be assured 
if a 6 is determined such that At ^ 6 for all likely nonstandard trajec- 
tories (At = t c - t n - At 0 ) . By defining AX such that AX = 0, t > t n + 5, 
the integral appearing on the right in equation (2,2.1) will always be 
zero and the second order term g 5 AX n At has been eliminated from the expres- 
sion for A0. 

At second stage ignition time, the only information available is 

AX°; “(t) will not be known for t > t 0 until a later time, and not 

completely known until cutoff. A value of X must be determined, however. 
With this in mind, AX will be defined as the sum of two functions. One 
function, AX 0 , will be determined at second stage ignition to meet the 

required end conditions if “(i) is zero throughout the second stage. 

The other function 5X will be determined as a function of Af/m as this 
function becomes known. Thus, AX is defined as 


AX = AX + SX. 

o 


( 2 . 2 . 2 ) 


Under the assumptions just described, AX 0 must satisfy the following 
equations . 


t n +S 


Ac = U-, AX + 
o 


/ 


t 

o 


V 5 


t n +5 


f x A> C q dt + / f 2 AX^dt + /f 3 AXq dt 


= 0, (2.2.3) 
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and 


t +8 t +5 t +8 

n n n 

dt + Jg 3 £Xg dt = 0. (2.2.4) 

t t t 

o o o 

The solution in the example given in Chapter I was obtained by 
assuming a constant coefficient system of differential equations for 
intervals of 40 seconds and the integral evaluated by summing over 
40— second intervals. This was done to simplify the problem to the point 
that the numerical results could be obtained by simple desk calculator 
operations. The numerical values employed in the following sections 
assumed a constant coefficient system over five-second intervals and 
evaluated the integral by summing over five-second intervals. This was 
done on an IBM 1620 digital computer and although differing only in the 
third significant figure in most cases, was considered necessary for the 
accuracy desired. 


A0 = Up AX + 


g x AX dt + 


§2 


From observation of the results of several different calculus of 
variations solutions, it was expected the following form would be ade- 
quate to represent AX q . 


AX 

o 


AC d - AC 2 t , 



(2.2.5) 


Since the standard trajectory is expected to always be included 
among the likely trajectories, 5 will always be less than or equal to 
zero. The value assigned to 8 in this report is -8.18 seconds. 

The requirements stated in equations (2.2,3) and (2.2.4) seem to 
imply that the simultaneous solution of two third-degree polynomials in 
two variables is needed. It should be remembered that if the linear 
terms are extremely large with respect to the higher order terms, the 
linear expression itself provides a good approximation to the solution. 
Such is the case in the problem at hand. The solution to the linear 
system can be used to approximate the second order term to effect an 
iteration on the solution. The linear system can be expressed as 
follows : 


BAC' = - U AX q , 


( 2 . 2 . 6 ) 


40 



where 


B 



flT 2 

glT 2 


dt, 


(2.2.7) 


and' 


U 



W(t n , 




AC' 


AC' 

O 

« 

A c 2 


( 2 . 2 . 8 ) 


(2.2.9) 


ACq and ACg, the solution to the linear system, are determined as follows 


AC' = - B" 1 U AX . 

o 


( 2 . 2 . 10 ) 


The second iteration gives 


AC" = 


AC" 

o 

AC'o 


from the following equation: 


BAC" 


+ U AX + Bo AC' = 0 


( 2 . 2 . 11 ) 
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where 


2AC f AC' 

O 2 


f 2 

f 2 T 2 

f 2 T 

§2 

§2T 2 

§2^ 


AC" = - B _1 U AX - B" 1 Bo AC' ' 


The third iteration including the third ordered term gives 


the solution to the following equation: 


BAC + U AX + B 2 AC" 2 + B 3 AC" 3 = 0, 


where AC" is obtained from equation (2.2.12) by replacing AC^ and AC' 0 
by ACq and AC", respectively, 
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AC" 3 


and 


3 

AC" 

o 

3 AC" 2 AC" 
o . 2 

3AC" AC" 2 
o 2 

_ Acf 


B 3 = 


t +& 
n 


t 

o 


f 3 

S3 


f 3 T 2 f 3 T 4 

g 3 T 2 g 3 T 4 


f 3 X 6 

g 3 T S 


dt. 


The solution for AC is 


1 (2.2.16) 


(2.2.17) 


AC = 


AC 

o 

ac 2 


B" 1 U AX q - B" 1 B 2 AC" 2 - B -1 B 3 AC" 3 . (2.2,18) 


The integrals and other numerical elements necessary to determine 
the matrices required of the above equations are listed in Appendix III. 
They have been obtained numerically under the assumption of a constant 
coefficient system over five-second intervals, and the integrals evaluated 
by summing over five-second intervals. These results yield the following 
matrices: 


Ui 


.349481 

1.259049 

.170980 

.491843 

U 2 


-.008175 

-.007129 

y 

-.005744 

-.008092 


(2.2.19) 


B 


-14.523714 

.570968 


-79.737064 

6.044026 


( 2 . 2 . 20 ) 
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-.092285 

-.328751 

-2.813252 

.003168 

.016477 

' .171907 


( 2 . 2 . 21 ) 


b 3 - 


.000737 

-.000029 


.004048 

-.000307 


.043109 .572059 

-.004554 -.075637 


( 2 . 2 . 22 ) 


Inverting the matrix in equation (2.2.20) gives 


- B - 1 


.143038 1.887071 

-.013512 -.343720 


(2.2.23) 


The result, together with equations (2.2.19), (2.2.21) and (2.2.23), 
gives 


- B -1 U = 


.034562 

-.001912 


.166639 .013617 .055082 

-.014562 -.000336 -.003864 


(2.2.24) 


-.007222 

-.015931 

-.078001 

.000158 

-.001221 

-.021075 


(2.2.25) 


and 


.000051 

0 

-.002427 

-.060906 

0 

.000051 

.000983 

.018268 


(2.2.26) 
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Choosing the following initial condition deviations, an example of 
the determination of AC will be followed. 


17.2860 


AX 

o 


5.5463 

-4.6494 


-58.34471 


(2.2.27) 


This, together with equations (2.2.10) and (2.2.24), gives 


AC' 


-1.7554 

.11319 


(2.2.28) 


The term AC' , defined in equation (2.2.12), is 


AC' 


2 


3.0814 
-.39739 . 


.012812 


This, with equations (2.2.14) and (2.2.25), gives 


-1.7554 


-.01692 


-1.7723 


4- 


= 


.11319 


.000702 


.11389 


2 3 

This result can be used to determine AC" and AC" as defined by 
equations (2.2.12) and (2.2.16), respectively, to give 


2 


AC" 


3.1410 

-.40369 

.012971 


and 


AC" 


-5.5668 

1.07320 

-.06897 

.001477 
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These results are used in equation (2.2.18), together with the 
matrices defined by equations (2.2.25) and (2.2.26), to give the 
following solution AC. 



-1.7554 


-.0173 


-.0002 

AC = 


+ 


+ 



.11319 


.00007 


.00001 


AC 


AC 

o 

AC 2 


-1.7727 
.11327 # 


(2.2.29) 


Section 3. Second Stage Perturbations 

Having obtained a function, AX 0 , which should meet the required end 
conditions for a standard second stage, it remains only to determine 6X 
such that the effect of Af/m on the end conditions is small. This 
requires further investigation of the expressions for Ar and AG which 
are obtained from equations (1.5.14) and (1.5.15). A restatement of 
these equations with the omission of the second order terms which have 
already been neglected or accounted for yields the following expressions. 


t 

n 


Ar - AX 4- 
x 0 



t 

o 


F AF dt. 


and 


(2.3.1) 


t 

n 


A0 - U 2 AX + 
^ o 


G AF dt. 


t 

o 


(2.3.2) 
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Substitution of U 2 for U x and G for F transforms the expression for Ar 
into the expression for A0. With this in mind, only the expression for 
Ar will be considered, and the corresponding result for A0 can be 
obtained by the substitution just described. From equation (2.2.2), 

AX = AX 0 + SX. Substituting this expression into AF in equation (2.3.1), 
and recalling that AF is defined as 


AF = 




1 + §rr ~ ) &X 2 


f/m 


1 + 0 ) 


Af/m 


(2.3.3) 


gives, after expanding and neglecting all previously unaccounted terms 
higher than second order, the following expression. 


Ar = 


U X £K + ■>- ! 

- 1 o / ± -o 


f i AX dt 4- / f P AX. 2 dt + 


p AAA. UU ~r I f 3 AX dt 

^ o o 


( £ i + 2f 2 ^ 0 > &x H f4 + ffe ^o) f + f2 5xS 


+ ■§+- — SX 
f/m m 


dt. 
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It will be recalled that the coefficients for £X 0 were determined 
to satisfy equation (2.2.3). This constraint required the following 
relationship to exist. 


Ui AX + 
- 1 o 


f, AX dt + 



fp ZNX 2 dt + 

o 


t 

o 


t 


o 



f 3 AX 2 dt = 0. 


Thus, the entire contribution to Ar stems from the last integral of the 
expression, namely, 


Ar 


rn r 


(f i + 2f, 


AX ) 

o 


SX + 





fp sx 2 


t 

o 


fi Af 
f/m m 


&X 


dt. 


(2.3.4) 


The evaluation of the above integral requires the knowledge of 
Af/m along the entire trajectory. This information is not available 
until after the cutoff conditions have been reached. However, in 
order that the integral be zero, it is sufficient that the integrand 
be zero everywhere. In determining an adequate guidance function, this 
restriction can be relaxed. It is sufficient that the integrand be 
sufficiently near zero that the integral is negligible . In order to 
. accomplish this, 5X must be a function of Af/m. The following form is 
chosen for 5X. 


SX 


= W-, — + w 2 

m ^ 



(2.3.5) 
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Substituting this expression into equation (2.3.4) and neglecting terms 
higher than second order gives 


n 


Ar = 


6 r i (t > IT + 5r ^ (t) (f 


dt, 


where 


Br 1 (t) = f jWi + f 4 + + 2f 2 AX, q W 


and 


Sr 2 (t) = f x W 2 + f 2 W x 2 + W x . 


Similarly, for A9, the following expression is obtained: 


n 


A0 = 


SerCt) f + Be a (t) 


2 


dt. 


where 


S9i(t) = 8l W x + g 4 + AK o + 2g 2 ZSC o W x 


and 


59 2 (t) = gl W 2 + g 2 W x 2 + jfe W x . 


(2.3.6) 


(2.3.7) 


(2.3.8) 


(2.3.9) 


.3.10) 


.3.11) 
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The term VJ ± will be considered a quadratic function of t and the 
coefficients determined to minimize the sum of squares of equations 
(2.3.7) and (2.3.10) simultaneously for a number of time points along 
the trajectory. W 2 will also be considered a quadratic function of t 
and the coefficients determined similarly using equations (2.3.8) and 
(2.3.11). For convenience, the following definitions will be used 
where the indicated summation is intended to include a convenient 
number of time points evenly distributed over the entire trajectory. 


1 


C(Br) = 


f , 2 



T 



(2.3.12) 


c(se) = 



(2.3.13) 


Ci(Sr) 



(2.3.14) 


1 


C3.(80) = 


y gig4 



(2.3.15) 


C L j(Sr ) = - 





(2.3.16) 


50 
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b = 


bi 

b 2 


( 2 . 3 . 21 ) 


b' = 


b 

b 

b 2 


( 2 . 3 . 22 ) 


C 4 (Sr) = - 


Ll 

f/m 


( 2 . 3 . 23 ) 


c 4 (se) = - 


V fiif 

L> f/m 


( 2 . 3 . 24 ) 


C 5 (Sr) - - / f 1^2 


( 2 . 3 . 25 ) 
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C 5 (68) - - ) § 1 S 2 


a = 


a o 

a i 

a 2 


(2.3.26) 


(2.3.27) 


With these definitions established, Wx and W 2 will be defined as 
follows and a least squares fit obtained. 


Wj. = b Q + b + b 2 x 2 . (2.3.28) 

W 2 = a Q + a x x + a 2 x 2 . (2.3.29) 


Minimizing the sum of squares of equation (2.3.7) over the 
specified number of time points yields the following constraint on the 
coefficients b 0 , b;L and b 2 : 


C(5r)b = C 1 (8r) + 


C 2 (5r) + 2b* (b) C 3 (&r) 


AC. 


The criterion for minimizing the sum of squares of 59 (t) defined 
in equation (2.3,10), however, is 


C(60)b = C 1 (50) + 


C 2 (S0) + 2b* (b) C 3 (60) 


AC. 
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In general, the least squares criterion for Sr^t) differs from 
that required for SGj^t). A single solution requires one criterion 
which is a function of the two separate criteria just derived. The 
single criterion used in this report will be a weighted sum of the two 
separate criteria just derived. To this end, the following matrices 
will be defined. 


D = C(6r) + Ct^CCSe). (2.3.30) 
Di = C x (5r) + a 2 C 1 (50) . (2.3.31) 
D 2 = C 2 (Sr) + a 2 C 2 (S0). (2.3.32) 
D 3 = C 3 (6r) + a 2 C 3 (60) . (2.3.33) 
D 4 = C 4 (8r) + a 2 C 4 (S0). (2.3.34) 
D 5 = C 5 (Sr) + a 2 C 5 (S0). (2.3.35) 


The symbol a is a weight used to represent the relative importance 
of angular error A0 to radial error, Ar. Combining the two separate 
criteria into one criterion in the manner described, yields the follow- 
ing constraint on b. 


Db = D x + 


D 2 + 2b* (b) D 3 


AC. 


The term in brackets represents a second order term and will be 
handled by obtaining an approximate solution b 1 by ignoring higher order 
terms and then determining b by using b f to approximate b in the expres- 
sion b*(b). With this in mind, the solution to the following equations 
defines b* * 


Db ? = D ! « 
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Then, 


b' = D" 1 D x . 


(2.3.36) 


With this first order approximation available, an iteration will be used 
to determine b from the equations 


Db 


Di + 


D 2 + 2b* (b* ) D 3 


AC, 


and 


b = D -1 D x + D" 1 


D h + 2b* (b' ) D 3 


AC, 


or 


b = b' + D" 1 


D 2 + 2b* (b 1 ) D 3 


AC. 


(2.3.37) 


Having determined the coefficients for W l5 it remains only to 
determine the coefficients for W 2 . Minimizing the sum of squares of 
6r 2 (t) gives 


C(8r)a = b* (b ) C 5 (Sr)b + C 4 (8r)b. 


Similarly for S0 2 (t), the following result is obtained: 
C(50)a = b* (b ) D 5 b + D 4 b. 

Combining these two expressions as before gives 
Da = b* (b ) D 5 b + D 4 b. 
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The contribution of W P (Af/m) 2 is small compared to that of W-, — . In 

m 

addition, b* is a good first order approximation to b, In view of these 
facts and to avoid unnecessarily complicated expressions introduced by 
relatively insignificant terms, b 1 will be used in the above system 
thus making "a" independent of AC. Under these conditions, the follow- 
ing equation will be used to determine "a". 


a 


= D 


-l 


- 

D 4 b' + b* (b 1 ) D 5 b ! . 


(2.3.38) 


Equations (2.3.36), (2.3.37) and (2.3.38) provide a means of 
determining the coefficients for W 4 and W p so as to minimize, as near 
as possible within the restrictions imposed, the sum of squares of 
Sr x ( t ) , 5r 2 (t), 56 x (t), and 50 p (t) an< 3 hence their effect on Ar and A0. 
An example will determine whether or not this has been accomplished 
sufficiently well for the result to be used as a guidance function. 

The values for the f’s and the g's were not available. The com- 
puter program employed yielded the f f s and g f s evaluated at five- 
second intervals and multiplied by At = 5 sec. Since every value used 
was multiplied by the same constant, the least squares solution is 
mathematically identical to that which would have been obtained if 
none of the values had been multiplied by At. Since the following 
summations were done on the desk calculator, only every fourth point 
was used and each summation involved 24 time points spaced 20 seconds 
apart beginning with t = 150 seconds. The following numerical values 
were obtained. 


25 


25 


25 


25 


25 



I 


fx 2 T = 



2_2 = 
1 T 




.662385 


1.112071 


2.782090 


8.040666 


25.186436 


25 


25 


25 


25 


25 


I 

I 

I 

I 

I 


\ 

gx 2 = .118972 x 10" 2 

gl 2 T = .427513 x 10- 2 

gl 2 T 2 = 1,664566 x 10“ 2 \ 

gl 2 r 3 = 6.737941 x 10- 2 

gl 2 T 4 = 27.921969 x 10‘ 2 

(2.3.39) 
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f if 4 = -3.417960 25 Y g x g 4 = - .26176 x 10 " 2 ^ 

fif 4 x = -4.104875 25 ^ g x g 4 T = - .55670 x 10 “ 2 )(2.3.40) 


V -1 


fif 4 x 2 = -7.961945 25 ) gl g 4 T 2 = -1.573233 x 10~7 




f T 2 

f/m 


.071016 25 


= .008462 x 10 -2 

f/m 


fl 2 T 


« .107982 25 


f/m 


fn 2 T 2 _ 
f/m 


fn 2 X 3 = 
f/m 


.656697 


f/m 


= .027354 x 10 " 2 

f/m 


= .244417 25 


2 2 i 

& 1- T - = .100396 x 10 " 2 } 

f/m 


(2.3.41) 


2^3 


25 ) = .390356 x 10 ~ 2 

f/m 


Z 2 4 

= 1.569404 x 10 -2 
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25 


25 


25 


25 


25 


25 


f 1 f 2 = .42248 x 10" 2 
^f!f 2 T = .59287 x 10" 2 
^ fif 2 T 2 = 1.27461 x 10“ 2 
£ fifsT 3 = 3.29295 x 10" 2 
7 fif 2 x 4 = 9.43480 x 10- 2 
y - 29.30378 x 10" 2 


25 y§i§2 


.040996 x 10- 


25 > gig 2 T = .102289 x 10 


-4 


25 ) gig 2 x^ = .316953 x 10 


-4 


25 ) gig 2 T 3 = 1.079346 x 10 


-4 


25 ) gig 2 ^ = 3.876986 x 10 


-4 i 


25 ) gig 2 ^ 5 = 14.529144 x 10 


— 4 


25 y fif2T 6 = 95.57426 x 10“ 2 25 ^gigsT 6 = 55.835969 x 10“ 4 

(2.3.42) 

To obtain one set of equations, instead of one for Ar and another 
for A9, it is necessary to choose a suitable value for CC. The value 
chosen for this report is a = 5. This is equivalent to saying that an 
error in angle A0 of 0.1 degree is five times as bad as an error in r, 
Ar, of 100 m. Or, equivalently, a value of Ar of 100m is weighted the 
same as A9 = .02 degrees. With a = 5, a 2 = 25. This, together with 
the sums evaluated in equations (2.3.39) through (2,3.42) gives the 
following matrices: 


25D 


.652128 

1.218949 

3.198232 


1,218949 

3.198232 

9.725151 


3.198232 

9.725151 

8.355253 


(2.3.43) 
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25D X = 


3.483400 

4.244050 

8.355253 


(2.3.44) 


25Dc. = 


.073132 -.269516 
.114820 -.754286 
-.269516 -2.334939 



- .432729 

- 1.353848 




- .618442 

• - 3.562787 



25D 3 = 

- 1.353848 

- 10.404047 

x 10" 2 



- 3.562787 

- 32.936066 




-10.404047 

-109.533252 




-.073132 

-.114820 

-.269516 


25D 4 = 

-.114820 

-.269516 

-.754286 



-.269516 

-.754286 

2.334939 



-.432729 

-.618442 

-1.353848 


-.618442 

-1.353848 

-3.562787 

25D 5 = 

-1.353848 

-3.562787 

-10.404047 


-3.562787 

-10.404047 

-32.936066 


-10.404047 

-32.936066 

-109.533252 


(2.3.45) 


(2.3.46) 


(2.3.47) 


x 10" 2 . (2.3.48) 
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From equation (2.3.43), the following inverse is obtained. 


25Dn 


= 25 D_1 = 


11.368704 

-11.105430 

2.227201 


•11.105430 

14.724394 

-3.347514 


2.227201 

•3.347514 

.821713 


. ( 2 .: 


Equations (2.3.36), (2.3.44), and (2.3.49) give the following 
solution for b' . 


b' = 


11.0786 

-4.1629 

.4168 


( 2 . 


Equations (2.3.45) and (2.3.49) yield 


D _1 D 2 = 


-.15656 

.02372 

.00002 


.11224 

-.29707 

.00607 


( 2 . 


Referring to equations (2.3.20), (2.3.46), (2.3.49) and (2.3. 
the following result is obtained. 


D' 1 b~(b') D 3 = 


11688 

.05583 

06228 

-.15880 

00861 

.03506 


( 2 . 


.49) 


3.50) 


3.51) 


50), 


3.52) 
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The results obtained from equations (2,3.50), (2.3,51), and 
(2,3.52) can be combined in equation (2,3,37) to give the following 
values for the coefficients of W x . 



b 

0 


11.0786 


-.39032 

.22390 

b = 

bi 

= 

-4.1629 

+ 

.14828 

-.61467 


_t*2 _ 


.4168 


-.01720 

.07619 


(2.3.53) 


From equations (2,3.47), (2.3,49) and (2.3.50), we get 


D 4 b = 


•1.7201 
.8397 
• .1083 


(2.3.54) 


Equations (2.3.20), (2.3,48), (2.3.49) and (2.3.50) give 


D" 1 b* (b ' ) D 4 b' 


-1.1589 
.7580 . 


- .11970 


(2.3.55) 


Using the expression obtained in equations (2.3,54) and (2.3.55) 
to evaluate equation (2.3.38) gives the following coefficients for W 2 . 



a 

0 


r 

1 


-2.8790 

a = 

a i 

a 2 

= D" 1 

D 4 b' + b*(b' ) D 5 b' 


1.5977 
- .2280 


(2.3.56) 
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The complete expression for X can now be stated. Equations (2.1.13), 
(2.2.2), (2.2.5), (2.3.5), (2.3.28), (2.3.29), (2.3.53), and (2.3.56) 
yield the following function. 


X = (58.2802 + AC 0 ) + 7.6268t + (.3574 + AC 2 )t 2 + (b 0 + b x x + b 2 T 2 ) 


+ (-2.8790 + 1 . 5977 t - .2280t 2 ) (Af/m) 2 , 


0 g t g 4.65685 


X = 58.2802 + 7.6268t + .3574t 2 , t > 4.65685 


(2.3.57) 


The terms AC Q and AC 2 are determined at second stage ignition from 
equations (2.2.10), (2.2.14), and (2.2.18) using the matrices evaluated 
in equations (2.2.24), (2.2.25) and (2.2.26) with the definitions of 
equations (2.2.9), (2.2.12) and (2.2.16). Also b G , b ± and b 2 are deter- 
mined at second stage ignition and are evaluated after AC c and AC 2 by 
equation (2.3.53). The symbol t is defined in equation (2.1.1), and 
Af/m is defined as follows. 




+ At 0 ) 


8.78065 

1.3751 - . 20888t ’ 


where Af/m is in m/sec 2 and X is in degrees. 


This guidance function has been derived for trajectories in the 
neighborhood of the standard trajectory employed. The example already 
used to illustrate the method of determining AX. Q employed the following 
initial condition deviations stated in equation (2.2.27) where Ax Q , Ay Q 
are in km and Ax 0 and Ay Q are in m/sec. 


AX 


Ax 

o 


17.2860 

Ay o 


5.5463 



-4.6494 



-58.34471 
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These values yielded the following values for AC 0 and AC 2 
(equation 2.2.29). 


AC = 


1 

o 

o 

<1 

1 


-1.7727 

1 

0) 

<] 

1 


.11327 


These values of AC 0 and AC 2 can be used in equation (2.3.53) to 
give the following values for b 0 , b x and b 2 . 



b 


11.7959 


o 



b = 

bi 

= 

-4.4954 


b 2 ^ 


.4559 


The function which would be used to determine X for the deviations in 
initial conditions used in this example would be 


X = 56.5075 + 7.6268T + .47067t 2 + (11.7959 - 4.4954t + ,4559t 2 ) — \ 

m | 

! 

+ (-2.8790 + 1.5977t - .2280t 2 ) (Af/m) 2 : 

/ 

0 S T S 4.65685 

r 

I 

X = 58.2802 + 7.6268t + .3574t 2 , t > 4.65685 ) 

(2.3.58) 


where 


Af = f 
m m 


(t + At c ) 


8.78065 

1.3751 - . 20888t * 
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Section 4, Implementation 


The employment of the guidance function just described would 
involve the precalculation of the following quantities which have been 
numerically evaluated for the mission under consideration. 


C Q = 58.2802 

C x = 7.6268 ) (2.4.1) 

C 2 = .3574 

J 


-B -1 U = 


.034562 

-.001912 


.166639 

-.014562 


.013617 

-.000336 


.055082 

-.003864 


(2.4.2) 


007222 

-.015931 

-.078001 

000158 

-.001221 

-.021075 


(2.4.3) 




.000051 

0 


0 

.000051 


-.002427 

.000983 


-.060906 

.018268 


(2.4.4) 


b' = D" 1 D 1 = 


11.0786 

-4.1629 

.4168 


^-l 


D 4 b' + b“ (b' ) D 5 b ' = 


-2.8790 
1.5977 . 
-.2280 


(2.4.5) 


(2.4.6) 
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With these quantities precomputed, they can be used to determine 
the following quantities as soon as AX 0 is determined. This is deter- 
mined by measuring x, y, x and y at ignition time of the second stage. 
From these measurements the elements of AX 0 are determined. 


X 


153.983 

o 



y o 


6435.878 

*o 


2818.329 



988.358 


(2.4.7) 


where x Q and y Q are measured in kilometers and x G and y Q are in m/sec. 
Then AC is computed in the following sequence of operations. 



.034562 

.166639 

.013617 

.055082 


AC' 

AC* = 





AX 0 = 

o 


-.001912 

-.014562 

-.000336 

-.003864 


AC 2 


(2.4.8) 



, 2 


r~ — 


AC' 


P *1 

-.007222 -.015931 -.078001 


O 

^ .2 


ACq 



2 AC' AC' 

= 


.000158 -.001221 -.021075 


O 

2 


ACg 



AC' 2 




(2.4.9) 
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and 


AC = AC' + 


-.007222 

.000158 


-.015931 

-.001221 


-.078001 

-.021075 


2 

AC" 

2AC'q AC'^ 

AC" 2 






ac'q 3 


.000051 

0 

-.002427 

-.060906 


3AC'q 2 AC" 


0 

.000051 

.000983 

.018268 


3AC" AC" 2 

O 2 





- 


.3 






AC" 



AC o 

AC 2 


(2.4.10) 

Having obtained the values of AC Q and AC 2 , the values of b Q , b± 
and b 2 can be determined. 


AC. (2.4.11) 


With b Q , b - l , b 2 , AC Q and AC 2 determined from the deviations at second 
stage ignition, AX 0 , the thrust angle X is determined hereafter by the 
foil owing f unc t ion . 


b 


11.0786 

0 



bi 

= 

-4.1629 

b 2 


.4168 


-.39032 

.22390 

.14828 

-.61467 

-.01720 

.07619 
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X = (58.2802 + AC q ) + 7.6268t + (.3574 + AC 2 ) t 2 + b Q ^ + b x T — N 


m 


+ b 2 T 2 — - 2.8790 (Af/m) 2 + 1.5977t (Af/m) 2 - .2280t 2 (Af/m) 2 , 
0 g t S 4.65685 


X = 58.2802 + 7.6268t + .3574x 2 , t > 4.65685 


(2.4.12) 


The term Af/m in the above expression is determined by the relationship 


Af , . f ( v 8.78065 

m ( ' t; m ‘ 1.3751 •- ,20888 t ’ 


(2.4.13) 


where 


T 


t* - t. 

1 

100 


(2.4.14) 


The term t^ is second stage ignition on the nonstandard trajectory, 
t 1 is time measured on the nonstandard trajectory; t' and t^ are in 
seconds, f/m is in m/sec 2 , and X is in degrees. 
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CHAPTER III 


RESULTS AND CONCLUSIONS 


Section 1, Adjustment to a Different Standard Trajectory 

In the event that performance of trajectories in the neighborhood 
of a different standard trajectory is to be investigated, it might be 
best' to duplicate the computations which have been described in Chapters 
I and II for the new standard trajectory if it differs greatly from the 
standard which was originally assumed. However, if it differs only 
slightly, there will be small deviations in the desired end conditions 
if the function derived for X is used to define X on this trajectory. 
These deviations will be represented by the following vector. 


AR = 


Ar 

A9 


( 3 . 1 . 1 ) 


To meet the desired end conditions, equation (2.2.15) will be modified 
to the following: 


AR = BAC + U AK q + B 2 AC" 2 + B 3 AC" 3 + AR = 0. (3.1.2) 

If the definition of AC 1 given by equation (2.2.6) is modified to the 
definition 


AC 1 = - B" 1 U AX - B" 1 AR, 
o 


(3.1.3) 


then all the subsequent computations can be performed exactly as before. 

To demonstrate the accuracy of the function which has been derived, 
an available computer program was employed which was designed for this 
purpose. The differential equations of this program described the 
motion in three dimensions and used a gravity field which differed from 
that assumed in this report. The program was extremely intricate and 
difficult to alter so that it would duplicate the simple two-dimensional 
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differential equations of motion which were assumed in deriving the 
guidance function presented in Chapter II. It was considered simpler 
to assume that this program represented a slightly different standard 
trajectory and alter the coefficients for X accordingly. The follow- 
ing values for Ar and A0 were obtained when X., defined by the coefficients 
in equation (2.1.13), was used. 


Ar = -1.8084 km 
A0 = .05955° 


and 


= 


-1.8084 

.05955 


(3.1.4) 


(3.1.5) 


Using -B -1 from equation (2.2.23), the following values are obtained. 


-B" 1 M = 


-.1464 

.0040 


Employing this result in equation (3.1.3), together with -B -1 U 
defined in equation (2.2.24), gives the following relationship for AC*. 


AC' 


.034562 

-.001912 


.166639 .013617 

-.014562 -.000336 


.055082 

-.003864 


^o 


+ 


-.1464 


.0040 


(3.1.6) 
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Section 2. Coefficient Computations 


Having determined the guidance function and adjusted it to the 
standard trajectory of an available computer program, it was then 
verified by using this guidance function for a number of different 
sets of initial conditions. These initial conditions resulted from 
various nonstandard first stages and were combined with several dif- 
ferent second stage perturbations. Table 3.1 lists these deviations 
in initial conditions along with the first stage deviation which caused 
them. Equation (3.1.6) was used to determine AC' instead of equation 
(2.4.8); AC was then determined from equations (2.4.9) and (2.4.10). 
These results along with b Q , b x and b 2 computed from equation (2.4.11) 
are listed in Table 3.2. 


TABLE 3.1 

DEVIATIONS IN INITIAL CONDITIONS 


Example 

No. 

First Stage 
Deviations 

Ax 

(km) 

Ay 

(km) 

Ax 

(m/ sec) 

Ay 

(m/sec) 

i 

None 

0 

0 

0 

0 

2 

+5000 lb 

-.93152 

-1.1623 

-30.2692 

-22.78749 

3 

-5000 lb 

.94444 

1.1824 

30.8335 

23.36393 

4 

Engine #2 
out at 100 sec. 

17.28599 

5.5463 

-4.6494 

-58.34471 

5 

Tail Wind 

.10183 

1.8180 

2.7670 

25.69483 

6 

Head Wind 

-.34235 

- .4055 

-5.3173 

-5.40216 

7 

Left Cross Wind 

-.10434 

- .1788 

-1.5250 

-2.33100 

8 

Right Cross Wind 

.09092 

.2155 

1.4240 

2.87015 

9 

-17= W 

4.76116 

1.1588 

48.5244 

-4.54684 

10 

+17= W 

-3.31692 

-.6959 

-39.7298 

3.74070 

11 

+17= F 

1.19032 

1.9245 

29.8569 

29.73483 

12 

-17= F 

-1.22157 

-1.9050 

-30.7597 

-29.22741 
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TABLE 3.2 


COEFFICIENTS DETERMINED FROM INITIAL CONDITIONS 


Example No. 

AC 0 

ac 2 

b o 

t>i 


i 

- .1464 

. 00400 

11.1381 

-4.1877 

.4197 

2 

-2.0636 

.12186 

11.9115 

-4.5439 

.4616 

3 

1.7735 

-.11491 

10.3634 

-3.8305 

.3777 

4 

-1.9224 

.11798 

11.8556 

-4.5206 

.4589 

5 

1.6001 

-.12227 

10.4293 

-3.8516 

.3801 

6 

- .5977 

.03329 

11.3205 

-4.2725 

.4297 

7 

- .3295 

.01634 

11.2122 

-4.2224 

.4238 

8 

.0702 

-.01088 

11.0504 

-4.1465 

.4148 

9 

.6194 

-.02062 

10.8343 

-4.0593 

.4047 

10 

- .7151 

.01946 

11.3631 

-4.2814 

.4306 

11 

2.2337 

-.15007 

10.1762 

-3.7408 

.3671 

12 

-2.5718 

.15881 

12.1178 

-4.6418 

.47 31 


Section 3. Results of Application 

Tables 3.3 and 3.4, respectively, show the deviations in Ar and A0 
obtained from these examples. Example No. 4 which resulted from an 
engine out at 100 seconds in the first stage had a rather large out-of- 
plane position and velocity deviation which was not assumed to exist 
in the equations for which the function was derived. In spite of this, 
the results reflect the type of accuracy that can be expected when mis- 
sion accomplishment is mathematically imposed as a criterion for deter- 
mining the coefficients of the guidance function. 
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Table 3.5 shows the additional fuel required beyond that required 
for the calculus of variations solution to the same set of conditions J 
No calculus of variations solution was available for the case in which 
F and W were each -5% throughout the second stage so that a fuel com- 
parison could not be made for this example. 


TABLE 3.3 
Ar (Meters) 


Perturbations 


2nd Stage 

1st Stage 

None 

-17c f, W 

+17c f, W 

-57c f, W 

None 

-i 

53 

-40 

23 

+5000 lb 

35 

95 

-27 

-157 

-5000 lb 

25 

55 

28 

151 

Engine out at 100 sec 

179 

175 

179 

-326 

Tail Wind 

* 

73 

7 

201 

Head Wind 

-2 

60 

-44 

-12 

Left Cross Wind 

-7 

54 

-43 

7 

Right Cross Wind 

* 

46 

-42 

32 

-1% W 

* 

61 

50 

* 

+17» 

* 

144 

-9 

* 

+1% f 

* 

57 

43 

* 

-17= f 

* 

92 

-21 

* 


*Due to underestimation of required computer time, these cases 
were not completed. It was felt that their inclusion would not 
significantly change the results. 
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TABLE 3.4 
A0 (Degrees) 


Perturbations 


2nd Stage 

1st Stage 

None 

-17= f, W 

+17= f, W 

-57= f, W 

None 

-.001 

.012 

-.027 

-.004 

+5000 lb 

.004 

.019 

-.026 

.020 

-5000 lb 

-.003 

.006 

-.023 

-.023 

Engine 2 out at 100 sec 

-.001 

.013 

-.029 

.008 

Tail Wind 

-.006 

.006 

-.030 

-.015 

Head Wind 

* 

.013 

-.026 

.000 

Left Cross Wind 

-.001 

.012 

-.027 

! -.002 

Right Cross Wind 

-.002 

.011 

-.027 

-.005 

-1% W 

* 

.009 

-.021 

* 

+1% W 

* 

.014 

-.032 

* 

+17= f 

* 

.005 

-.023 

'k 

-17= f 

* ! 

.021 

-.025 

* 
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TABLE 3.5 


Fuel Loss (lbs) 


Perturbations 


2nd Stage 
1st Stage ^ 

None 

-17, f, W 

+17o f, W 

None 

5 

5 

5 

+5000 lb 

7 

7 

6 

-5000 lb 

4 

4 

4 

Engine 2 out at 100 sec 

2 

1 

2 

Tail Wind 

5 

4 

5 

Head Wind 

JU 

5 

5 

Left Cross Wind 

5 

5 

5 

Right Cross Wind 

5 

5 

5 

-1% W 

* 

5 

5 

+1% W 

* 

7 

7 

+1% f 


4 

4 

-17o f 

* 

7 

6 


Note: Calculus of variations solutions for the -5% f, W case 

were not available to compare fuel consumption. 
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It is interesting that the only requirements which would keep fuel 
consumption small was the fact that X was fitted to the value of X obtained 
from the calculus of variations solution for the standard trajectory and 
the AX 0 was chosen of a form which would fit well the deviations • 
encountered by calculus of variations solutions with nonstandard initial 
conditions. Aside from these two requirements, no further attempt was 
made to minimize fuel. Nevertheless, Table 3.5 indicates that this 
appears to be sufficient for negligible fuel loss. 

To illustrate the difference in some of these nonstandard trajec- 
tories, it should be pointed out that example No. 4 which resulted from 
a first stage engine out at 100 seconds with second stage perturbations 
of -57 0 thrust deviation and simultaneously -5% flow rate deviation 
required 30 seconds additional burning time and displaced the cutoff 
point 136 km. The radius error was -326 m and cutoff angle error +.008 
degrees both of which appear quite acceptable. 


Section 4. Further Applications , 

The guidance function derived and applied in this report was pre- 
sented primarily to illustrate the effectiveness of obtaining an explicit 
solution to the linearized differential equations of motion. It should 
be emphasized that to linearize the differential equations, it is not 
necessary to linearize every parameter. It is sufficient that the 
state variables be linear leaving the forcing functions to whatever 
form they may have. This same procedure can be applied to other dif- 
ferential equations including calculus of variations equations or equa- 
tions of motion in a different coordinate system. The linearized equa- 
tions are always sufficiently accurate in some neighborhood of the 
standard solution. Whether this neighborhood is sufficiently large to 
include all expected deviations cannot be stated in general but must be 
specifically investigated for the particular differential equations under 
consideration. The differential equations of motion are readily adaptable 
to such analysis. 

In addition to the parameters included in the analysis of Chapter I, 
any of a number of other parameters could have been included and their % 
effect on cutoff deviations determined. 


Section 5. Conclusions 


The results of this investigation illustrate the accuracy of the 
solution to the linearized differential equations of motion as well as 
the value and means of obtaining the solution in explicit form. In 
addition to the insight it provides, it is a powerful tool for mathe- 
matically imposing the mission criteria in the determination of the 
coefficients of a guidance function. Further applications of the 
explicit solution to a linearized set of differential equations to the 
analysis of a nonlinear system are restricted only by the differential 
equations themselves and the imagination and ingenuity of the analyst. 
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APPENDIX I 


U(t n , t) 

(t n = 620.68 sec.) 


U(t n , 


1 

0 

0 

0 


0 0 0 
1 0 0 
0 10 
0 0 1 


U(t n , 600) 


U(t n , 560) 


1 

0 

.02608 

0 

0 

1 

0 

.02608 

-.01886 

.02833 

1 

0 

.02833 

.04812 

0 

1 


.99917 

.00104 

.06068 

.00000 

.00104 

1.00200 

.00000 

.06068 

-.05890 

.07863 

.99925 

.00113 

.07863 • 

.14472 

.00113 

1.00192 



U(t n , 520) 


U (t n , 480) 


U(t n , 440) 


U(t n , 400) 


.99647 

.00366 

-.10332 

.12186 

.99177 
.00760 
-.14971 
_ .16106 

.98488 

.01239 

-.19814 

.19521 

.97571 

.01767 

-.24801 

.22453 


.00366 

1.00808 

.12191 

.24516 

.00760 

1.01846 

.16133 

.34912 

.01240 

1.03332 

.19605 

.45658 

.01771 

1.05283 

.22654 

.56761 


.10065 

.00004 

.99689 

.00428 

.14051 

.00019 

.99276 

.00915 

.18018 

.00049 

.98677 

.01577 

.21958 

.00099 

.97884 

.02340 


.00004 

.10076 

.00428 

1.00771 

.00019 
.14108 
.00916 
1 . 01752 - 

.00049 

.18182 

.01561 

1.03148_ 

.00099 
.22315 
.02345 
1 .04974 
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.96418 

.02326 

.25861 

.00170 

U(t n , 

360) = 

.02315 

-.29892 

1.07716 

.25337 

.00170 

.96892 

.26526 

.03251 



.24933 

.68251 

.03238 

1.07244 



.95028 

.02880 

.29718 

.00263 

U(t n , 

320) = 

.02855 

-.35047 

1.10648 

.27702 

.00263 

.95696 

.30835 

.04264 



_ .26978 

.80158 

.04235 

1.09974 



.93395 

.03412 

.33519 

.00378 


.03362 

1.14106 

.00377 

.35261 

U(t n , 280) = 

-.40248 

.29797 

.94294 

.05372 


.28603 

.92547 

.05314 

1.13180 



.91513 

.03904 

.37255 

.00514 

U(t n , 240) = 

.03813 

1.18117 

.00511 

.39825 

-.45485 

.31671 

.92684 

.06564 


.29820 

1.05475 

.06458 

1.16882 


U(t n , 200) 


U(t n , 160) 


U(t n , t Q ) 


.89380 

.04182 

-.50739 

.30618 

.86987 
. 04446 
-.56010 
.30993 

.86115 

.04503 

-.57740 

.31006 


.04336 

1.22718 

.33355 

1.19028 

.04693 

1.27957 

.34893 

1.33309 

.04792 

1.29901 

.35385 

1.38267 


.40916 

.00664 

.90865 

.07651 

.44491 

.00831 

.88835 

.08876 

.45637 

.00890 

.88097 

.09284 


.00670 

.44550 

.07831 

1.21101 

.00843 

.49459 

.09165 

1.25862 

.00905 

.51145 

.09625 

1.27619 
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APPENDIX II 


5(t j) 



Uii 

U 12 (10“ 2 ) 

u X3 (io- 4 ) 

fi 14 

160 

.8296 

-1.259 

-.4212 

12.69 

200 

.8888 

-1.548 

-.4512 

14.64 

240 

.7704 

-1.549 

-.3911 

13.67 

280 

.6480 

-1.534 

-.3290 

12.58 

320 

.5228 

-1.501 

-.2654 

11.37 

36Q 

.3968 

-1.447 

-.2015 

10.06 

400 

.2768 

-1.367 

-.1385 

8.65 

440 

.1556 

-1.258 

-.0790 

7.17 

480 

.0508 

-1.106 

-.0258 

5.62 

520 

-.0308 

- .904 

.0156 

4.02 

560 

-.0748 

- .630 

.0380 

2.41 

600 

-.0533 

- .254 

.0271 

.81 

600* 

-.0426 

- .203 

.0216 

.81 


/V a 

The elements of U(t.) were determined by multiplying the elements of 
U(t n , tj) H(tj) by At j . Except for t = 160 and t = 600, Atj = 40 sec 
was used. For t = 160, Atj = 33.19 sec. and for t = 600, At^. = 40.68. 
For t = 600*, Atj = 32.5 sec was used for the first three columns since 
their multiplier is a power of AX which, for the guidance function, was 
defined to be zero for the last 8.18 seconds. These values should be 
used instead of those at 600 in evaluating the integral for the guidance 
function. 
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t . 

,1 

U 2i 

U 22 (10- 2 ) 

u 23 (io- 4 ) 

U 24 

160 

-1.585 

-.895 

.805 

9.014 

200 

-1.900 

-.897 

.965 

8.475 

240 

-1.874 

-.762 

.952 

6.730 

280 

-1.834 

-.630 

.931 

5.168 

320 

-1.775 

-.500 

.901 

3.791 

360 

-1.696 

-.374 

.861 

2.603 

400 

-1.590 

-.254 

.807 

1.610 

440 

-1.452 

-.144 

.737 

.821 

480 

-1.272 

-.048 

. 646 

.242 

520 

-1.037 

.026 

.526 

-.117 

560 

- .723 

.065 

.367 

-.249 

600 

- .296 

.046 

.150 

-.147 

600* 

- .237 

.037 

.120 

-.147 


'See Footnote on Page 83. 
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tj 

U31 

U 32 (10" 2 ) 

u 33 (io- 4 ) 

U 34 

160 

1.414 

- 2.676 

-.718 

26.960 

200 

1.701 

-3.563 

-.864 

33.676 

240 

1.666 

-3.951 

-.846 

34.879 

280 

1.598 

-4.390 

-.811 

36.005 

320 

1.486 

-4.888 

-.754 

37.036 

360 

1.319 

- 5.457 

-.670 

37.944 

400 

1.081 

-6.117 

-.549 

38.708 

440 

.749 

-6.898 

-.380 

39.338 

480 

.290 

-7.816 

-.147 

39.690 

520 

-.347 

-8.949 

.176 

39.855 

560 

-1.242 

-10.383 

.631 

39.754 

600 

-2.583 

-12.286 

1.311 

39.364 

600* 

-2.064 

-9.816 

1.048 

39.364 


See Footnote on Page 83. 
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t . 

,1 

041 

u 42 (io- 2 ) 

u 43 (io- 4 ) 

u 44 

160 

-3.904 

-2.468 

1.982 

24.868 

200 

-5.035 

- 2.657 

2.555 

25.112 

240 

-5.396 

-2.442 

2.739 

21.561 

280 

-5.804 

-2.209 

2.947 

18.121 

320 

-6.271 

-1.950 

3.184 

14.778 

360 

-6.815 

-1.656 

3.460 

11.514 

400 

-7.456 

-1.314 

3.785 

8.312 

440 

-8.226 

- .906 

4.176 

5.167 

480 

-9.173 

- .405 

4.657 

2.055 

520 

-10.369 

.229 

5.264 

- 1.018 

560 

-11.933 

1.063 

6.058 

- 4.071 

600 

-14.318 

2.216 

7.270 

- 7.100 

600* 

-11.439 

1.770 

5.808 

- 7.100 


See Footnote on Page 83. 
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APPENDIX III 


NUMERICAL RESULTS FROM 5 SECOND INTERVALS 


Limits 

Integrand 

t 

n 

I 

t o 

t n +b 

f 

fc n _ . 

fc n +5 

I 

_ 

fi 

-14.541427 

.017713 

-14.523714 

flT 

-29.213315 

.082930 

-29.130385 

flT 2 

-80.125334 

.388272 

-79.737062 

f x T 3 

-251.962610 

1.817852 

-250.144758 

flT 4 

-857.620677 

8.511001 

-849.109676 

flT 5 

-3076.064786 

39.847657 

-3036.217129 

f lT 6 

-11,454.32276 

186.56275 

-11,267.76001 

f 2 

-.092309 

.000024 

-.092285 

f 2 T 

-.142486 

.000113 

-.142373 

f 2 x 2 

-.329282 

.000531 

-.328751 

f 2 X 3 

-.914112 

.002487 

-.911625 

f 2 x 4 

-2.824898 

.011646 

-2.813252 

f 2 x 5 

-9.372182 

.054524 

-9.317658 

f 2 X 6 

-32.711274 

.255278 

-32.455996 

f 3 

.000738 

-.000001 

.000737 

f 3 T 

.001483 

-.000004 

.001479 

f 3 x 2 

.004068 

-.000020 

.004048 

f 3 X 3 

.012792 

-.000092 

.012700 

f 3 x 4 

.043541 

-.000432 

.043109 

f 3 x 5 

.156170 

-.002023 

.154147 

f 3 x 6 

.581531 

-.009472 

.572059 
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Limits 

Integrand 

t 

n 

f 

t 

0 

t +5 

n 

/ 

t 

n 

t +6 
n 

r 

j 

t 

n 

§i 

.594373 

-.023405 

.570968 

81T 

1.824835 

-.109580 

1.715254 

SlT 2 

6.557068 

-.513042 

6.044026 

gl^ 3 

25.223896 

-2.402011 

22.821885 

gl^ 4 

100.949256 

-11.245976 

89.703280 

gl^ 5 

414.750164 

-52.652539 

362.097625 

gl^ 6 

1736.333387 

-246.513923 

1489.819464 

gs 

.003170 

-.000002 

.003168 

g2^ 

.006143 

-.000012 

.006131 

gs^ 2 

.016530 

-.000055 

.016475 

ga^ 3 

.051330 

-.000256 

.051074 

g2^ 4 

.173072 

-.001197 

.171875 

g2^ 5 

.615988 

-.005604 

.610384 

g2T 6 

2.278515 

-.026239 

2.252276 

S 3 

-.000030 

.000001 

-.000029 

g 3 ? 

-.000093 

.000006 

-.000087 

g 3T 2 

-.000333 

.000026 

-.000307 

g3T 3 

-.001281 

.000122 

-.001159 

g 3^ 4 

-.005125 

.000571 

-.004554 

g3-l' 5 

-.021057 

.002673 

-.018384 

g3T G 

-.088153 

.012516 

-.075637 







349481 

1.259049 

.170980 

.491843 

.008175 

-.007129 

-.005744 

-.008092 
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